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Abstract: We study the time evolution of unstable dS p x S q configurations with flux in 
theories of gravity with a cosmological constant. For certain values of the flux, we identify a 
stable configuration to which these unstable solutions flow. For other values of the flux the 
sphere wants to decompactify, regardless of the sign of the initial perturbation. 
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1. Introduction 

It has been known since the work of Freund and Rubin [1] that when an antisymmetric tensor 
field of dimension <7 — 1 is added to gravity (Einstein-Hilbert action), space-times of dimension 
p + q can naturally compactify to product spaces of the form dS p x S q . 

Bousso, DeWolfe and Myers [2] derived analogous results when a positive cosmological 
constant is added to the gravitational action, and studied the stability of these compactifica- 
tions under gravitational perturbations. Their analysis showed the stability to depend on the 
relative value of the flux compared to the cosmological constant as well as on the dimension 
of the internal space. 

The unstable modes, at these compactifications, were identified in [2], generalizing the 
work of [3]. In this paper, we derive their equations of motion and study their evolution. In 
doing so we will make the assumption that these remain the only unstable modes throughout 
the evolution. A similar analysis, in the absence of flux, was recently done by Contaldi, 
Kofrnan and Peloso [4], 

This paper answers the questions about the fate of these configurations that was posed 
in the original paper of Bousso, DeWolfe and Myers. It is unclear to what extent the lessons 
learned here translate into more realistic models, e.g. [5, 6]. In fact, Giddings and Myers [7] 
have already studied these types of models and argued that positive vacuum energy together 
with extra dimensions render the four-dimensional space unstable toward decompatification 
of the extra dimensions. Our study is purely classical and does not incorporate the effect of 
thermal fluctuations or tunneling considered in [7]. 
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After this work was completed, we learned that the time evolution of these configurations 
had already been studied in the past [8, 9]. Our results present a more complete analysis 
of all the initial conditions and agree with them when these conditions overlap. The goal of 
these earlier papers was to use the unstable mode as the inflaton mode. Other references that 
partially overlap with this work are [10],[11],[12]. 

In section 2 we present a summary of the solution described in [2] to establish notation 
and state the problem. In section 3, we display a more general solution of the equations 
of motion that reduces to the solutions found in [2] and hence can interpolate between an 
unstable initial configuration and a stable final state. In section 4, we will present the time 
evolution that results from the numerical analysis. 

2. Product Spacetimes with Flux 

We will consider solutions of the following action; 

S = IJ - 2A - ^ F Pl ... Pq F Pl " Pq ) (2.1) 

in units where Mp in p + q dimensions has been set equal to 1. The cosmological constant 
A will be assumed positive. Fp l p q is a totally antisymmetric tensor of rank q. Throughout 
this paper we adopt the notation that upper case latin indices run from 1 to p + q, greek 
indices run from 1 to p and lower case latin indices run from 1 to q. The equations of motion 
are 


Gmn = 


1 


2( g -l)! 


Fmp 1 ...p„- 1 F 1 


Pi...P q -i 

N 


- 9mn (A + -^F Pl ... Pq F Pl - Pq ) (2.2) 


d M (V=9F Mp i- Pq -')=0 (2.3) 

These equations admit a solution that is a product of two spaces, a Lorentzian and a Rie- 
mannian ( S q ): 


1 

7 =9 


ds 2 = —dt 1 + a{t) z d p + R^dQ 


(2.4) 


F q = (P 


W 0 I 59 j) V 0 I 59 = i?Q Qq (2-5) 

To be consistent with the notation of [2] we will parametrize the solutions in terms of 


T = — 
4A 


Though T has been often referred to, in the literature, as the flux, it is clear from (2.5) that 
the actual flux, with the normalization given above, is also a function of the compactification 
radius Rq. With a normalization that makes the flux radius independent, F q = Rq ^ volgg, 
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there are actually two static solutions, instead of one, for some values of c. The equations of 
motion for the Lorentzian scale factor, a(t) are: 


(P~ 1)(P~2) a 2 (p- l)(p~2) k = q{q- 1) 1 

2 a 2 2 a 2 1 1 2 R 2 

( 2 . 6 ) 

k = 1,0, —1 corresponds to closed, flat or open spaces. The effective cosmological constant of 
the p dimensional Lorentzian space, A p . is given by: 

(P - !) A p = ~ , 2 r f_ 2 (1 - (? - W (2-7) 

P I Z 

Since A has been chosen positive, the effective cosmological constant will be positive as long 
as 


The radius of S q is given by: 


0 < 7 < - 

q- 1 


( 2 . 8 ) 


q -i 
R 2 


2A 


(l + (p-l):F) 


(2.9) 


p + q ~ 2 

The following table summarizes the regions of stability as a function of the free dimensionless 
parameter 7 and the dimension of the Riemannian space q as found in [2], 



dS p x S q 
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q 
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Not reflected in the previous table is the fact that the number of unstable modes increases 
as a function of q. For q = 2, and q = 3 there is only one unstable mode. It corresponds to 
the radius of S q . For q > 4 other modes became tachyonic as well. 

3. Evolution of Unstable Spacetimes 

There is a broader set of solutions to (2.2),(2.3), that do not correspond to a product space- 
time. They are of the form: 


ds 2 = — dt 2 + a(t) 2 d p x x + R(t) 2 dVt q 


(3.1) 
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vols"? = R{t) q Qq 


F = _ i _ 

q W(t) 


V 0 I 59 


(3.2) 


In this case, the normalization of the differential form F q is fixed by the equations of motion. 
The corresponding equations of motion are: 


qjq - 1 ) 
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A ( X A 

(3.3) 


Q f 2 . 

where ^. To lighten these already complicated equations we have made the assumption 
that fc = 0, the generalization to k = ±1 is straightforward. These equations can be reduced 
into an equation for the breathing mode of the sphere, often called the radion: R(t) 


h _q + p-2R 2 ^p-l h 
p-2 R p — 2 ' 

where the potential 


lq(p + q ~ 2) R 2 2 P ~ 2 
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The evolution of a(t ) is given by: 
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(3.4) 


(3.5) 


(3.6) 


The initial conditions, i.e. the sign of - at /?, = 0, and i? equal the radius of the unstable 
point, decide which branch of the equations (3.4), (3.6) to consider. 

As long as the flux is nonvanishing and relatively small, 0 < f3 s , where 


_ (q - l) 2q 1 fp + q - 2\ q 

P~ 1 V 2( 1 J 


(3.7) 


the radion potential will always have two stationary solutions, one stable and one unstable. 
There is a solution to the equations of motion that corresponds to the field R(t) sitting on 
the maximum. In this case, the geometry factorizes and corresponds to a dS p x S q space. 
There are also solutions that correspond to R(t) sitting in the minimum. In these cases, the 
corresponding geometry also factorizes into the product of the spaces L p x S q . The Lorentzian 
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Figure 1: Potential for the field R(t) as a function of the flux, f3, when p = 4, and q = 2. 


signature space will be a (A)dS or flat Minkowski depending on the value of (3. In particular, 
for (3 m < (3 < p s , where 


(3m 


(q — l) 2q 1 
29 


(3.8) 


the Lorentzian signature space will also be a de Sitter space. 

We are interested in studying the evolution of the unstable solutions away from the 
maximum. The equations (3.4),(3.6) will capture such an evolution provided that R(t ) is 
indeed the only mode excited along the path away from the maximum. The analysis of [2] 
proves that this is indeed the only excited mode around the extrema of the potential, provided 
that q = 2,3, for q > 4 there will be other unstable modes. In this work we will restrict our 
attention to the first two cases, and leave the analysis of q > 4 to future work. 


4. Results 

In this section we report on the numerical solution of the equations (3.4), (3.6). It is useful 
to distinguish between two cases: 

• When the minimum of the potential corresponds to another de Sitter solution (/3 m < 

(3 < (3 S ) 

• When the minimum of the potential corresponds to an Anti de Sitter solution (0 < (3 < 

Pm) 
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Figure 2: Value of ^ at the stationary points as a function of the initial condition /?. The upper 
branch corresponds to the unstable point and the lower branch corresponds to the stable point. (Case 

P = 4, q = 2.) 
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Figure 3: Numerical evolution for p = 4, q = 2, f3 = 0.3. The arrows indicate the direction of time 
evolution. 

4.1 P m < P < Ps 

Let us choose the initial condition such that a/a > 0, namely, we start with an expanding 
de Sitter space. The numerical analysis shows that the evolution from the unstable dS p x S q 
in this case leads either to a decompactification of the compact dimensions, or to another, 
stable dS p x S q solution, see Figure 3. 

In the case of decompactification, the expansion rate of the radius of the sphere S q 
asymptotes to: 


R I 2A 
R ~ y (D-l)(D-2) 
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Figure 4: Flux (3 as a function of parameter T . when p = 4 and q = 2. 

where D = p + q- It is interesting to note that this expansion rate matches the asymptotic 
effective Hubble constant for the de Sitter space, i.e. R/R = a/a as R, a —> oo. This agrees 
with the result found in [4], because this fast evolution erases all the information about the 
initial flux. 

The other direction of evolution from the unstable dS p x S q , with these initial conditions, 
leads to a stable dS p x S q solution. The two solutions are characterized by the same value of 
/3. Namely, for a static solution we have: 

P/A q = RR 2q 

2 _ (g-!)(£>- 2) 

2 A i ~( P - i)r 

For a particular value of (3 we have two static solutions related by 

tj f i + (p- 1)^> y 

Ti v 1 + ip - l )Fi ) 

The stable solution will always have a smaller radius R than the initial, unstable one. The 
two are related by: 



(4.3) 


(4.2) 



1 + (P ~ l )Fi t>2 

l + ( P -l)Tf * 


(4.5) 


The effective cosmological constant for the de Sitter is given by (2.7), and it therefore changes: 


f 2 _ 1 + (f - q)Ff ( «y 2 
\o/ / 1 + (1 - q)Ti \a) l 


(4.6) 
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Figure 5: The evolution of effective Hubble rate for p = 4, q = 2. (The colors match the numerical 
evolution in Figure 3) 


In particular, for p = 4 and q = 2 we have: 


Tf 9 Ti 

(4.7) 

p2 _ ^ p2 _ o'TT p2 

K f ~ K i — oJ-iKi 

(4.8) 

fd \ 2 (9^-1) / a \ 2 
\a) f 95i(l-7i) Wi 

(4.9) 

. A4 i — rpjK 

A4 f = " A 

1 6A 4i - A 

(4.10) 


Notice, that the final effective cosmological constant can be made very small compared 
with the initial. An analogous transition has been found in the context of braneworlds [13]. 
The radius of the internal sphere will remain, however, within the same order of magnitude 
of the initial one. Let us compare the entropies of the initial and the final solution: 


S{dS 4 x S 2 ) f _ A{dSi) f x V (S' 2 )/ _ (1-^X1 +3 T/) 2 1 -ft 

S{dS 4 xS 2 )i A(dSi)i x V {S 2 )i (1- ^)(l + 3^) 1 l) 9ft - 1 


as long as, | which corresponds to our initial unstable static solution. 

Now let’s look at the situation where the initial conditions are chosen such that the de 
Sitter component of the space is initially contracting, a/a < 0 (See Figure 6). In this case, 
the numerical evolution shows that the internal dimensions necessarily decompactify, while 
the initial de Sitter dimensions keep contracting. This result may appear surprising, at first, 
when the initial velocity at the unstable point is taken to be negative, i.e. R < 0. These 
initial conditions choose the negative branch in (3.6) and hence the positive branch in (3.4). 
The friction force in (3.4) is then big enough to overcome the force generated by the potential. 
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Figure 7: Numerical evolution for p = 4, q = 2, (3 = 0.15 or equivalently Anti de Sitter minimum. 

The crunching solution should (3.1) asymptotically become a Kasner type solution. We 
did not use this ansatz in the numerical analysis that led to Figure 6. 


4.2 0 < (3 < p m 

When the flux j3 is very small, only one of the static solutions corresponds to a de Sitter 
compactification, the other solution corresponds to an anti de Sitter space. In this case the 
numerical solution in both branches (either initially contracting (Figure 8), or expanding 
(Figure 7) de Sitter phase) leads to a decompactification of the inner space, while the initially 
de Sitter components either expand or contract, depending on the initial conditions and the 
initial direction of the numerical evolution. This is a sensible result. The classical trajectories 
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Figure 8: Numerical evolution for p = 4, q = 2, (3 = 0.15 or equivalently Anti de Sitter minimum. 

studied here cannot connect spaces with different spacial curvatures. Given the ansatz for 
the metric (3.1) the initial de Sitter configuration can be either flat or closed while the final 
anti de Sitter configuration has to be open. 

5. Conclusions 

In this paper we have studied the evolution of gravitationally unstable dS p x S q compactifi- 
cations. We have identified an additional static stable solution that only exists for a certain 
interval of fluxes. Under certain conditions the evolution leads to this static solution. The 
other possible final states correspond to a decompactifying internal sphere, while the de Sitter 
space may either expand or contract. 
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